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Abstract 



We present the currently most accurate evaluation of the W boson mass, Mw, 
in the Minimal Supersymmetric Standard Model (MSSM). The full complex phase 
dependence at the one-loop level, all available MSSM two-loop corrections as well as the 
full Standard Model result have been included. We analyse the impact of the different 
,^ , sectors of the MSSM at the one-loop level with a particular emphasis on the effect of 

the complex phases. We discuss the prediction for Mw based on all known higher-order 
contributions in representative MSSM scenarios. Furthermore we obtain an estimate of 
the remaining theoretical uncertainty from unknown higher-order corrections. 
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1 Introduction 

The relation between the massive gauge-boson masses, M^ and Mz, in terms of the Fermi 
constant, G^, and the fine structure constant, a, is of central importance for testing the elec- 
troweak theory. It is usually employed for predicting M^ in the model under consideration. 
This prediction can then be compared with the corresponding experimental value. The cur- 
rent experimental accuracy for Mw, obtained at LEP and the Tevatron, is 6Mw = 30 MeV 
(0.04%) [1,2]. This experimental resolution provides a high sensitivity to quantum effects 
involving the whole structure of a given model. The M^-Mz interdependence is therefore 
an important tool for discriminating between the Standard Model (SM) and any extension 
or alternative of it, and for deriving indirect constraints on unknown parameters such as the 
masses of the SM Higgs boson or supersymmetric particles, see Ref. [3] for a recent review. 
Within the SM the confrontation of the theory prediction and experimental result for the 
W boson mass supplemented by the other precision observables yields an indirect constraint 
on the Higgs-boson mass, Mh, of Mh = 89ttl GeV with Mh < 175 GeV at the 95% C.L. [2]. 
Within the Minimal Supersymmetric Standard Model (MSSM) [4] the electroweak precision 
observables exhibit a certain preference for a relatively low scale of supersymmetric particles, 
see e.g. Refs. [5,6]. 

The experimental precision on Mw will further improve within the next years. The 
Tevatron data will reduce the experimental error to about 6Mw = 20 MeV [7], while at the 
LHC an accuracy of about 6Mw = 15 MeV [8] is expected. At the GigaZ option of a linear 
e"^e~ collider, a precision of 6Mw = 7 MeV can be achieved [9, 10]. 

A precise theoretical prediction for Mw in terms of the model parameters is of utmost 
importance for present and future electroweak precision tests. Within the SM, the com- 
plete one-loop [11] and two-loop [12-16] results are known as well as leading higher-order 
contributions [17-21]. 

The theoretical evaluation of Mw within the MSSM is not as advanced as in the SM. So 
far, the one- loop contributions have been evaluated completely [22-26], restricting however 
to the special case of vanishing complex phases (contributions to the p parameter with non- 
vanishing complex phases in the scalar top and bottom mass matrices have been considered 
in Ref. [27]). At the two-loop level, the leading 0{aas) corrections [28,29] and, most recently, 
the leading electroweak corrections of 0{a^), 0{atab), 0{al)to Ap have been obtained [30, 
31]. Going beyond the minimal SUSY model and allowing for non-minimal flavor violation 
the leading one-loop contributions are known [32]. 

In order to fully exploit the experimental precision for testing supersymmetry (SUSY) and 
deriving constraints on the supersymmetric parameters,^ it is desirable to have a prediction 
of Mw in the MSSM at the same level of accuracy as in the SM. As a step into this direction, 
we perform in this paper a complete one-loop calculation of all contributions to Mw in the 
MSSM with complex parameters (cMSSM), taking into account for the flrst time the full phase 
dependence and imposing no restrictions on the various soft SUSY-breaking parameters. We 
combine this result with the full set of higher-order contributions in the SM and with all 
available corrections in the MSSM. In this way we obtain the currently most complete result 
for Mw in the MSSM. A public computer code based on our result for Mw is in preparation. 



"'^ A precise prediction for Mw in the MSSM is also needed as a part of the "SPA Convention and Project" , 
see Ref. [33]. 



We analyse the numerical results for M^ for various scenarios in the unconstrained MSSM 
and for SPS benchmark scenarios [34]. The dependence of the result on the complex phases 
of the soft SUSY-breaking parameters is investigated. We estimate the remaining theoretical 
uncertainties from unknown higher-order corrections. 

The outline of the paper is as follows. In Sect. 2 the basic relations needed for the 
prediction of M^/ are given and our conventions and notations for the different SUSY sectors 
are defined. In Sect. 3 the complete one-loop result Mw including the full phase dependence 
of the complex MSSM parameters is obtained. The incorporation of all available higher-order 
corrections in the SM and the MSSM is described. The numerical analysis is presented in 
Sect. 4, where we also estimate the remaining theoretical uncertainties in the prediction for 
M\Y. We conclude with Sect. 5. 

2 Prediction for Mw — basic entries 

Muons decay to almost 100% into eUf^Ue [35]. Historically, this decay process was first de- 
scribed within Fermi's effective theory. The muon decay rate is related to Fermi's constant, 
Gfj,, by the defining equation 



with F[x) = 1 — 8x — 12x^lnx + 8x^ — x^. By convention, the QED corrections in the 
effective theory, Aqed, are included in eq. (1) as well as the (numerically insignificant) term 
3m^/(5M^) arising from the tree-level W propagator. The precise measurement of the muon 
lifetime and the equivalently precise calculation of Aqed [36,37] thus provide the accurate 
value 

G^ = (1.16637 ± 0.00001 x 10"^) GeV^l (2) 

In the SM and in the MSSM, G^ is determined as a function of the basic model parameters. 
The corresponding relation can be written as follows, 

G e^ 

^ = —, ^^ (1 + Ar). (3) 

The quantity Ar summarizes the non-QED quantum corrections, since QED quantum effects 
are already included in the definition of G^ according to eq. (1), which makes the evaluation 
of Ar insensitive to infrared divergences. Ar depends on all the model parameters, which 
enter through the virtual states of all particles in loop diagrams, 

Ar = Ar{Mw,Mz,mt,a,as,...,X) (4) 

with 

X = Mhsm (SM), 

X = Mh, Mh, Ma, Mh±, tan/5, Mj, Af, m^o.±, . . . (MSSM), 



and is hence a model-specific quantity. For each specification of the parameter set X, the 
basic relation (3) is fulfilled only by one value of M]^. This value can be considered the 
model-specific prediction of the W mass, Mw[X), providing a sensitive precision observable, 
with different results in the SM and in the MSSM. Comparing the prediction of M\y with 
the experimental measurement leads to stringent tests of these models. Within the SM 
one can put bounds on the Higgs boson mass, which has so far not been directly measured 
experimentally. In the MSSM, Ar and Mw are sensitive to all the free parameters of the 
model, such as SUSY particle masses, mixing angles and couplings. The SUSY entries are 
briefly described below, thereby introducing also conventions and notations for the subsequent 
discussion. 

Sfermions: 

The mass matrix for the two sfermions of a given flavour, in the fi, Jr basis, is given by 

/M| + mf mfX}\ 



with 



Ml = M| + M| cos 2/3 (J/ -Q^s^ 



F 
r2 



M^ = M|, + Mi COS 2/3 g^s4 (6) 

Xf = Af — fi* {cot (3 ,tanP}, 

where {cot P , tan/3} applies for up- and down-type sfermions, respectively, and s^ = sin^ ^w = 
1 — My^r/M"^. In the Higgs and scalar fermion sector of the cMSSM, Nf + 1 phases are present, 
one for each Af and one for /i, i.e. Nf + 1 new parameters appear. As an abbreviation, 

<pAf = arg(A/) , (7) 

will be used. As an independent parameter one can trade (pAf for (pXf = arg(X/). The 
sfermion mass eigenstates are obtained by the transformation 

h ) '\f. 

with a unitary matrix U?. The mass eigenvalues are given by 

1 



and are independent of the phase of Xf. 

Higgs bosons: 

Contrary to the SM, in the MSSM two Higgs doublets are required. At the tree level, the 
Higgs sector can be described with the help of two independent parameters, usually chosen 
as the ratio of the two vacuum expectation values, tan/3 = V2lv\, and M4, the mass of the 
CP-odd A boson. Diagonalisation of the bilinear part of the Higgs potential, i.e. the Higgs 



mass matrices, is performed via rotations by angles a for the CP-even part and /3 for the 
CP-odd and the charged part. The angle a is thereby determined through 



tan 2a = tan 2/5 



Ml +M2 



vr 



M\ 



MV 



< a < 

2 



(10) 



One obtains five physical states, the neutral CP-even Higgs bosons /i, iJ, the CP-odd Higgs 
boson A, and the charged Higgs bosons H^ . Furthermore there are three unphysical Gold- 
stone boson states, G'^.G^. At lowest order, the Higgs boson masses can be expressed in 
terms of M^, Mp^, and Ma-, 



M\ 



h,H 



Mi± 



MX + MIt \ {Ml + M|)2 - 4M2M| cos2 2/3 



Ml + M? 



w 



'li: 



Higher-order corrections, especially to M/j, however, are in general quite large. Therefore we 
use the results for the Higgs boson masses as obtained from the code FeynHiggs2.2 [38- 
40]; also the angle a is replaced by the effective mixing angle aefr in the improved Born 
approximation [41,42]. 



Charginos and neutralinos: 

The physical masses of the charginos are determined by the matrix 



X 



Ms y2sin/3MvK 

\/2 cos (3 Myy H 



(12) 



which contains the soft breaking term M2 and the Higgsino mass term /i, both of which may 
have complex values in the cMSSM. Their complex phases are denoted by 



^2 



arg (M2) and 



arg (/i) . 



(13) 



The physical masses are denoted as m-± and are obtained by applying the diagonalisation 
matrices Uj^± and V^±, 



u;±xvt± 



^x± 







^x^ 



(14) 



The situation is similar for the neutralino masses, which can be calculated from the mass 
matrix (sw = sin^^, c^ = cos9^) 



/Ml 

M2 

— Mz s^cos/3 Mz Cw cos /3 

y Mz Sw sin (3 Mz c^ sin (3 



-Mz s„ cos l3 Mz s„ sin l3\ 

Mz Cw cos (3 Mz Cw sin (3 

-/i 

-/i / 



(15) 



This symmetric matrix contains the additional complex soft-breaking parameter Mi, where 
the complex phase of Mi is given by 



mi 



arg (Ml). 



(16) 



The physical masses are denoted as m^o and are obtained in a diagonahsation proce- 
dure using the matrix N^o. 



NtoYNl 





Xi 
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t 





X2 








X" 








m^o 







V 








m^J 



(17) 



At the two-loop level also the gluino enters the calculation of M^/. In our calculation below 
we will incorporate the full phase dependence of the complex parameters at the one-loop 
level, while we neglect the explicit dependence on the complex phases beyond the one-loop 
order. Accordingly, we take the soft SUSY-breaking parameter associated with the gluino 
mass, M3 = rrig, which enters only at two-loop order, to be real. 

3 Calculation of Ar 

Our aim is to obtain a maximally precise and general prediction for Mw in the MSSM. So 
far the one-loop result has been known only for the case of real SUSY parameters [24,25]. In 
this section, we evaluate the complete one-loop result in the cMSSM with general, complex 
parameters and describe the incorporation of higher-order terms. 

3.1 Complete one-loop result in the complex MSSM 

Evaluation of the full one-loop results requires renormalisation of the tree-level Lagrangian. 
This introduces a set of one-loop counter terms, which contribute to the muon decay ampli- 
tude, in addition to the genuine one-loop graphs. At one-loop order, Ar can be expressed in 
terms of the W boson self-energy, vertex corrections ("vertex"), box diagrams ("box"), and 
counterterms for charge, mass, and field renormalisation as follows, 

Ar(°) = 



S^^(O) 




+ 2^- 
e 


si 


(5 Ml 
\Ml 




+\5Zl + 


\5Zl^ 


+ \5Z^L 


+ \5zl^ + 


(vertex) 



(box), (18) 

where Sr(g^) denotes the transverse part of a vector boson self-energy. The leading contribu- 
tions to Ar*^"^ arise from the renormalisation of the electric charge and the weak mixing angle 
(the last two terms in the first line of eq. (18)). The former receives large fermionic contribu- 
tions from the shift in the fine structure constant due to light fermions, Aa oc log(m//M^), 
a pure SM contribution. The latter involves the leading universal corrections induced by the 
mass splitting between fields in an isospin doublet [43], 

E?(0)_3-(0) 

In the SM Ap reduces to the well-known quadratic term in the top-quark mass if the masses of 
the light fermions are neglected. In the MSSM Ap receives additional sfermion contributions. 



in particular from the squarks of the third generation. The one-loop result for Ar expressed 
in terms of Aa and Ap reads 

Ar(") = Aa- %Ap + Ar^^l, (20) 

where Arrem summarizes the remainder terms in eq. (18). 

Throughout our calculation the on-shell renormalisation scheme is applied. In this scheme 
renormalisation conditions are imposed such that the particle masses are the poles of the 
propagators and the fields are renormalised by requiring unity residues of the poles. These 
conditions ensure that the Mw^Mz interdependence given by eq. (3) is a relation between 
the physical masses of the two gauge bosons. Applying the renormalisation conditions, the 
counterterms in eq. (18) can be expressed in terms of self-energies, and eq. (18) turns into 

w) 



^'" = ^^^ + (vertex) + (box) ^^j^ 



+ 



-^Re 



fc2=o Cw M| si 



Er(M|) E^'^iM^) 



Ml M^ 



-S2(0)-S^(0)-Sl^(0)-S^-(0), (21) 

where Si(g^) denotes the left-handed part of a fermion self-energy. The electron and muon 
masses are neglected in the fermion field renormalisation constants, which is possible since 
the only mass-singular virtual photon contribution is already contained in Aqed of eq. (1) 
and is not part of Ar. 

At the one-loop level one can divide the diagrams contributing to Ar into four classes: 

(i) SM-like contributions of quark and lepton loops in the gauge-boson self-energies, schemat- 
ically depicted in Fig. 1; 

(ii) SUSY contributions of squark and slepton loops in the gauge-boson self-energies, de- 
picted in Fig. 2; 

(iii) contributions from the Higgs and gauge boson sector, which contain besides self-energies 
(Figs. 3 and 4) also vertex and box graphs (Fig. 5). 

(iv) SUSY contributions involving neutralinos and charginos in self-energies, vertex graphs 
and box diagrams, see Figs. 6 and 7. 

The calculation of the one-loop diagrams was performed in the dimensional regularisa- 
tion [44] as well as the dimensional reduction scheme [45] . The analytical result for Ar turned 
out to be independent of the choice of scheme. 

Technically, the calculation was done with support of the Mathematica packages FeynArts 
[46,47], OneCalc [48], and FormCalc [49]. A peculiarity among the box diagrams is the graph 
with a virtual photon, see diagram (4) in Fig. 5. Since QED corrections are accounted for 
in the Fermi model definition of G^, according to eq. (1), the corresponding diagram with 
a point-like W^-propagator 1/M^ has to be subtracted. This is the same procedure as in 
the SM and leaves an IR-finite expression for Ar (see Ref. [11] for the original analysis, a 




(1) 

Figure 1: Generic lepton and quark one- loop diagram contributing to Ar via gauge-boson 
self-energies; Vi, V2 = 7, Z, W^, /i, /2 = i^, I, u, d. 
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Figure 2: Generic squark and slepton one-loop diagrams which contribute to Ar via gauge- 
boson self-energies; /, /i, /2 = u, I, u, d. 
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Figure 3: Generic gauge-boson contributions to one-loop gauge-boson and fermion self- 
energies entering Ar (the same diagrams as in the SM). The labels / and u in the fermion 
self-energy diagram stand for electron, muon and the corresponding neutrinos. 
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Figure 4: Generic contributions of MSSM Higgs bosons and Goldstone bosons to one-loop 
gauge-boson self-energies; s, Si, S2 = h, H, A, H"^, G^, G^. 
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Figure 5: Generic gauge-boson contributions to one-loop vertex and box diagrams entering 
Ar (the same diagrams as in the SM). The labels / and u in the vertex diagrams stand for 
electron, muon and the corresponding neutrinos. 
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Figure 6: Generic neutralino/chargino contributions to gauge-boson (1) and fermion (2) self- 
energy diagrams. The labels /, v in the fermion self-energy diagram stand for electron, muon 
and their corresponding neutrinos, and the labels /, v indicate their respective superpartners. 
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Figure 7: Generic neutralino and chargino contributions to one-loop vertex and box diagrams 
entering Ar. The labels /, v stand for electron, muon and their corresponding neutrinos, and 
the labels /, z> indicate their respective superpartners. 



detailed discussion of this point can be found in Refs. [12,13]). After this subtraction, all 
external momenta and all lepton masses can be neglected in the vertex and box diagrams, so 
that this set of Feynman diagrams reduces to one-loop vacuum integrals. 

In order to obtain the relevant contributions to Ar, the Born- level amplitude needs to 
be factored out. This is straightforward for the SM vertex and box graphs and also for the 
SUSY vertex corrections. Concerning the box contributions, diagrams like graph (3) of Fig. 5 
directly yield a Born-like structure (after simplifying the Dirac chain using, for instance, the 
Chisholm identity, 7^7j.7p = -ie^yp^'^'''^^ + g^y'^p - g^p'^^ + gvpl^i), i-e. 



M 



Born-likc box 



M 



Born 



Arbo 



(22) 



Here 

Mboth = {Uu^Y i^-Uf,){UelpUJ_V^J ■ ^ ^ (23) 

is the tree-level matrix element in the limit where the momentum exchange is neglected. 
Box graphs involving supersymmetric particles, on the other hand, yield a different spinor 
structure. The SUSY box graphs shown in Fig. 7, diagrams (3) and (4), can schematically 
be written as 

A^SUSYboxl = {UeYuJ^U^){u^^-fpUJ^V^J ■b'f\ (24) 

-A/isusYbox2 = {u^,uj^Ui,){u^^uj+Ve) ■ b'f\ (25) 

where eq. (24) corresponds to diagram (3) and eq. (25) to diagram (4). The SUSY box 
contributions to Ar can be extracted by applying the Fierz identities 

{UeY^-U^^){u^^lpUJ-V^J = -{u^^Y ^-U^,){UelpUJ-V^J (26) 



where eq. (27) can be further manipulated using charge conjugation transformations. This 
yields 

— Oi , ArsusY box2 — — ; 

lira Aira 



-^^SUSYboxl — Ol , ^^SUSY box2 — — ; O2 • \^^ 



3.2 Incorporation of higher-order contributions 

In order to make a reliable prediction for M^^^ in the MSSM, the incorporation of contributions 
beyond one-loop order is indispensable. We now combine the one-loop result described in 
the previous section with all known SM and MSSM higher-order contributions. In this way 
we obtain the currently most accurate prediction for My\/ in the MSSM. 

3.2.1 Combining SM and MSSM contributions 

As mentioned before, the theoretical evaluation of M^^ (or Ar) in the SM is significantly more 
advanced than in the MSSM. In order to obtain a most accurate prediction for M^y (via Ar) 
within the MSSM it is therefore useful to take all known SM corrections into account. This 
can be done by writing the MSSM prediction for Ar as 

where Ar^^ is the prediction in the SM and Ar'^^^^"^^ is the difference between the MSSM 
and the SM prediction. 

In order to obtain Ar'^^^^ according to eq. (29) we evaluate Ar'^^^^"^^ at the level of 
precision of the known MSSM corrections, while for Ar^^ we use the currently most advanced 
result in the SM including all known higher-order corrections. As a consequence, Ar^^ takes 
into account higher-order contributions which are only known for SM particles in the loop but 
not for their superpartners (e.g. two-loop electroweak corrections beyond the leading Yukawa 
contributions and three-loop corrections of (^(aa^)). 



It is obvious that the incorporation of all known SM contributions according to eq. (29) 
is advantageous in the decoupling limit, where all superpartners are heavy and the Higgs 
sector becomes SM-like. In this case the second term in eq. (29) goes to zero, so that the 
MSSM result approaches the SM result with M^sm = M^, where M^ denotes the mass of the 
lightest CV-even Higgs boson in the MSSM. For lower values of the scale of supersymmetry 
the contribution from supersymmetric particles in the loop can be of comparable size as the 
known SM corrections. In view of the experimental bounds on the masses of the supersym- 
metric particles (and the fact that supersymmetry has to be broken), however, a complete 
cancellation between the SM and supersymmetric contributions is not expected. Therefore 
it seems appropriate to apply eq. (29) also in this case (see also the discussion in Ref. [30]). 

3.2.2 SM contributions 

As mentioned above, within the SM the complete two-loop result has been obtained for 
Mw [12-16]. Besides the one- loop part of Ar [11] it consists of the fermionic electroweak 
two-loop contributions [12], the purely bosonic two-loop contributions [13] and the QCD 
corrections of 0{aas) [14,15]. Higher-order QCD corrections are known at 0{aal) [17,18]. 
Leading electroweak contributions of order 0{G'ias'mf) and 0{G^m^) that enter via the 
quantity Ap have been calculated in Ref. [20]. Furthermore, purely fermionic three- and four- 
loop contributions were obtained in Ref. [50], but turned out to be numerically very small 
due to accidental cancellations. The class of four-loop contributions obtained in Ref. [51] 
give rise to a numerically negligible effect. 

All numerically relevant contributions were combined in Ref. [16], and a compact ex- 
pression for the total SM result for M^ was presented. This compact expression approx- 
imates the full SM-result for Mw to better than 0.5 MeV for Higgs masses ranging from 
10 GeV < Mh < I TeV, with the other parameters (m*, Aal^{Mz), as(M^), Mz) varied 
within 2a around their central experimental values. The contributions entering the result 
given in Ref. [16] can be written as 

where we have suppressed the index "SM" on the right-hand side. 

3.2.3 MSSM two-loop contributions 

The leading SUSY QCD corrections of 0{aas) entering via the quantity Ap arise from 
diagrams as shown in Fig. 8, involving gluon and gluino exchange in (s)top-(s)bottom loops. 
These contributions were evaluated in Ref. [28]. We have incorporated this result into the 
term Ar^^^SM-SM j^ ^^ (29). 

Besides the 0{aas) contributions, recently also the leading electroweak two-loop cor- 
rections of 0{af), 0{al) and 0{atab) to Ap have become available [30]. These two-loop 
Yukawa coupling contributions are due to MSSM Higgs and Higgsino exchange in the (s)top- 
(s)bottom-loops, see Fig. 9. In Ref. [30] the dependence of the 0{affj) corrections on the 
lightest MSSM Higgs boson mass, M^, has been analysed. Formally, at this order the ap- 
proximation Mh = would have to be employed. However, it has been shown in Ref. [30] 
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t,b_ 



t,b 




(g) 



i~9) 



Figure 8: Sample diagrams for the SUSY 0{aas) contributions to Ap: {g) squark loop with 
gluon exchange, (g) (s)quark loop with gluino exchange. 

how a non-vanishing MSSM Higgs boson mass can be consistently taken into account, in- 
cluding higher-order corrections. Correspondingly we use the result of Ref. [30] for arbitrary 
Mfi and employ the code FeynHiggs2 . 2 [38-40] for the evaluation of the MSSM Higgs sector 
parameters. 




(q) 



(q) 



(H) 



Figure 9: Sample diagrams for the three classes of MSSM 0{a'^), 0{al), 0{atab) contribu- 
tions to Ap: (g) quark loop with Higgs exchange, (g) squark loop with Higgs exchange, (H) 
quark/squark loop with Higgsino exchange, (f denotes Higgs and Goldstone boson exchange. 

The irreducible supersymmetric two-loop contributions to Ar discussed above need to be 
supplemented by the leading reducible two-loop corrections. The latter can be obtained by 
expanding the resummation formula [53] 



1 + Ar 



1 



;i-A«)(l + ^Ap)-Ar, 



(31) 



up to the two-loop order. At this order eq. (31) correctly incorporates terms of the type 
(Aa)^, (Ap)^, AaAp and AaArj-em- In this way we account for the leading terms of order 
0(iVrQ;^), where Nf is the number of fermions. 

The final step is the inclusion of the complex MSSM parameters into the two-loop re- 
sults. So far all two-loop results have been obtained for real input parameters. We therefore 
approximate the two-loop result for a complex phase 0, Mp|^"(0), by a simple interpolation 
based on the full phase dependence at the one-loop level and the known two-loop results for 



11 



full 



real parameters, M^ 



(0), M, 



full, 
W ' 



vr 



M^\<P) = M'^{<P) + [M^\0)-M^{0)] 



1 



COSi 



+ 



rlL/ 



[M'^\n) - Mi^in 



2 
1 — cos < 



(32) 



Here M^{(f)) denotes the one-loop result, for which the full phase dependence is known. The 
factors involving cos0 ensure a smooth interpolation such that the known results M^^^(O), 
it) are recovered for vanishing complex phase. In Sect. 4.3 we estimate the uncertainty 



M, 



full, 
W 



in the Mw prediction due to this approximate inclusion of the complex phase dependence at 
the two-loop level. 

3.3 Practical determination of Ar 

As can be seen from eqs. (3), (4) M^ is directly related to Ar, which however depends on M^y 
itself. In practical calculations we therefore find M]y, using eq. (3), in an iterative procedure. 
According to our strategy outlined above we obtain M]y from Ar^^ and Ar^^^^~^'^. In order 
to make contact with the known SM result we use the compact expression for the total SM 
result for M^y as given in Ref. [16]. This requires some care in the iterative evaluation of M^. 
Since the compact expression for the SM result gives M^, its inversion yields Ar^^(M^) 
rather than Ar^^ as function of the MSSM value of Mw, which should be inserted in eqs. (3), 
(29). 

The desired expression for Ar^^{Mw) is approximately given by 



Ar^'^'iMw) 



Ar*^"^ + Ar*^"""^ + Ar^"°^^ + Ar*^"^^™'^ + A 



+ Ar^*^^""™'^ + Ar*^"^^™'^ 



(M, 



+ 



Ar 



fcrm, sublcad ' bos 



W 



{Mf^^ 



W )) 



(«2) 
ferm, red 



(33) 



where Ar 



(a2) 
fcmi, red 



are reducible two-loop contributions arising from eq. (31). All terms of the 



first two lines of the right-hand side of eq. (33) are known analytically. The numerically small 
contribution in the third line can therefore be obtained as 



Ar 



ferm, sublcad 



-F Ar 



bos 



{M^^) 



ArS^^(MS^^; 



Ar*^"^ -|- Ar'-""''-' + Ar*^""^^ 



_|_ ^j^iGlmf) _^ ^, 



,("') 



fcmi, red 
3„,6x 



_|_ /\.i:{Gf,(^sTnj) _|_ ^^(G;im^) 



(M^'] 



(34) 



from Ar^^{Ml^) [16]. Here A^ 



ferm, sublead 



denotes the subleading fermionic two-loop terms 



beyond the leading m'f term and the reducible terms. In this way the correct M\y dependence 
in Ar is only neglected in the subleading electroweak two-loop contributions Ar^^^^ subiead + 
Ar|j"g , which are numerically small. For Ar*^"^ we use our recalculation, while for Ar*^""''^ 
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^^.{aa^)^ ^^(G^mj^ ^^(G^Qsnij) ^^^^ /\fiG-^"h) we use compact expressions from 
Refs. [18,20,30,52]. Our final result for Ar^^^SM^^^) ^^^^^g 

where Ar^^{Mw) is given in eq. (33) and Ar^^^^~^^{Miv) is the difference between the full 
MSSM one-loop result as described in Sect. 3.1 and the SM one- loop result supplemented 
by the higher-order supersymmetric contributions specified in Sect. 3.2.3. Inserting this ex- 
pression into eq. (3) one can now calculate M]y using a standard iteration which is rapidly 
convergent. 

4 Numerical analysis 

In the following subsections we present our numerical results. First the impact of the one- 
loop contributions from the different sectors of the MSSM is systematically analysed, and the 
dependence of the result on the different masses and complex phases is studied in detail. As 
a second step we take into account all higher-order corrections and discuss the full prediction 
for the W boson mass in the MSSM for a choice of sample scenarios. 

For the numerical analysis the analytical results for Ar and Mw, which were calculated as 
described above, were implemented into a Fortran program. Though built up from scratch, 
for the calculation of the MSSM particle spectrum our code partially relies on routines which 
are part of the FormCalc package [49]. The Higgs sector parameters are obtained from 
the program FeynHiggs2.2 [38-40]. This Fortran program for the calculation of precision 
observables within the MSSM will be made publicly available [54]. 

If not stated otherwise, in the numerical analysis below for simplicity we choose all soft 
SUSY-breaking parameters in the diagonal entries of the sfermion mass matrices to be equal, 

Mj = Mp = Mp, = . . . , (36) 

see eq. (5). In the neutralino sector the GUT relation 

Ml = ^'-fM2 (37) 

(for real values) is often used to reduce the number of free MSSM parameters. We have 
kept Ml as a free parameter in our analytical calculations, but will use the GUT relation to 
specify Mi for our numerical analysis if not stated otherwise. 
We have fixed the SM input parameters as ^ 

G^ = 1.16637 X 10~^ Mz = 91.1875 GeV, as{Mz) = 0.117, 

a = 1/137.03599911, Aa[f2i = 0.02761 [55], Aaiep = 0.031498 [56], (38) 

mt = 172.5 GeV m^ = 4.7 GeV, rrir = rric = . . . = 



2t 



For the results shown in Sect. 4.1 we use mt ~ Yll.l ^ 2.9 GeV [57]. For the comparison of theoretical 
predictions with experimental data in Sect. 4.2 we use the most up to date value mt — 172.5 ± 2.3 GeV, 
which has recently become available [58] . 
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The complex phases appearing in the cMSSM are experimentally constrained by their contri- 
bution to electric dipole moments of heavy quarks [59] , of the electron and the neutron (see 
Refs. [60,61] and references therein), and of deuterium [62]. While SM contributions enter 
only at the three-loop level, due to its complex phases the cMSSM can contribute already at 
one-loop order. Large phases in the first two generations of (s)fermions can only be accomo- 
dated if these generations are assumed to be very heavy [63] or large cancellations occur [64] , 
see however the discussion in Ref. [65]. Accordingly (using the convention that 0m2 = 0), 
in particular the phase 0^ is tightly constrained [66], while the bounds on the phases of the 
third generation trilinear couplings are much weaker. 

4.1 Analysis of parameter and phase dependence 

We begin by studying the impact of the one-loop contributions to Ar from the various MSSM 
sectors, i.e. the sfermion sector, the chargino and neutralino sector, as well as the gauge boson 
and Higgs sector. In order to be able to analyse the different sectors separately, we do not 
solve eq. (3) using our complete result for Ar, but we rather investigate the mass shift 6M\y 
arising from changing Ar by the amount Ar^^^"^, 

SMw = - ^f /- , Ar^^^^. (39) 

Here Ar^^^^ represents the one-loop contribution from the supersymmetric particles of the 
considered sector of the MSSM. We fix the (in principle arbitrary) reference value for Mw in 
eq. (39) to be M^^ = 80.425 GeV. Our full result for My/, which is determined from eq. (3) 
in an iterative procedure, is of course independent of this reference value. 

4.1.1 Sfermion sector dependence 

We first investigate the influence of sfermion one-loop contributions, which enter via the 
gauge-boson self-energy diagrams depicted in Fig. 2. The selectron and smuon contributions 
to the electron and muon field renormalisations shown in Fig. 6 and to the vertex and box 
diagrams shown in Fig. 7 will be discussed as part of the chargino and neutralino contributions 
in Sect. 4.1.2. 

The leading one- loop SUSY contributions to Ar arise from the t/h doublet. Since the 
mass of the partner fermion appears in the sfermion mass matrices, see eq. (5), a significant 
splitting between the diagonal entries can be induced in the stop sector. The off-diagonal 
elements in the stop sector and for large tan/? also in the sbottom sector can furthermore 
give rise to a large mixing between the two states of one flavour. 

The complex parameters in the t/h sector are /i. At and Ai,. Neither the /i nor the A 
parameters appear explicitly in the couplings of the diagrams of Fig. 2. They only enter via 
the absolute values and phases of Xt^b-, the off-diagonal entries of the squark mixing matrices. 
We have checked at the analytical level that the phases (pXt t drop out entirely in the full 
one- loop calculation of Ar and have no influence on Myy. Hence, the phases and absolute 
values of /x. At and Af, enter the sfermion-loop contributions (at one-loop order) only via 

|Xip = |AP + lAicot/5 p -2|A| ■ |At|cot/? cos(0At +0m), (40) 
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In particular, the phases (pAt b ^^^ ^(j. only enter in the combinations {(j)^^ ^ + 0^) and only 
via modifications of the squark masses and mixing angles. 

The phase dependence is illustrated in Figs. 10 and 11, where the squark loop contri- 
butions to 6Mw (evaluated from eq. (39)) are shown as function of the phase combination 
{(pA + 0^) with (pAt = (pAb- Since the phases enter only via \Xtfi\, their infiuence is most 
significant if all terms in eqs. (40), (41) are of a similar magnitude. This is the case if tan/? is 
rather small and |/i| and \At^b\ are of the same order. Such a situation is displayed in Fig. 10 
and the left panel of Fig. 11, where tan/5 = 5 and \At\ = \Ai,\ = 2M? has been chosen. Fig. 10 
shows the effect on 5Mv/ from varying the phase (0a + 0^) for a fixed value of |/i| = 900 GeV 
and Mr = 500, 600, 1000 GeV, while in Fig. 11 the squark sector contributions to 5M\y are 
shown as contour lines in the plane of (0^ + 0^) and |/x|. In the scenario with tan/3 = 5 
(Fig. 10 and left panel of Fig. 11) the variation of the complex phase (0A + 0/i) can amount to 
a shift in the W boson mass of more than 20 MeV. The most pronounced phase dependence 
is obtained for the largest sfermion mixing, i.e. the smallest value of M? and the largest value 

of |/i|. 

The right panel of Fig. 11 shows a scenario where tan/5 is rather large, tan/5 = 30. As 
a consequence, \Xt\ ~ \At\ and \Xi,\ ^ |/itan/5|, so that the absolute values of Xt and X^ 
depend only very weakly on the complex phases. The plot clearly displays the resulting much 
weaker phase dependence compared to the scenario in the left panel of Fig. 11. The variation 
of the complex phase gives rise only to shifts in M^y of less than 0.5 MeV, while changing |/i| 
between 100 and 500 GeV leads to a shift in M^ of about 2 MeV. 

In the following plots we discuss the dependence of the sfermion loop contributions to 
SM\Y on the common sfermion mass M?. In Fig. 12 we show the squark contributions as 
function of M? for various values of {(pAt + 0/i) and for (0A(, + 0^) = 0. The intermediate 
value tan/5 = 10 is chosen, and \At^b\ = 350 GeV, \fi\ = 300 GeV. In agreement with the 
discussion above, the phase dependence is relatively small. It leads to a shift of about 5 MeV 
in Mw for low values of M? ^ 250 GeV, where the stop mixing is large. The total squark 
contributions can shift the prediction of M^y by up to 30 MeV for small M?. For large M? the 
squark contributions show the expected decoupling behaviour. However, even for sfermion 
masses as large as M? = 1000 GeV the shift in Mw is still about 4 MeV, i.e. half the size of 
the anticipated GigaZ accuracy. 

In Fig. 13 we show the squark and slepton contributions for the same parameters as before 
except that {(pAt + (pp.) = and {(pA^ + <Pfi) is varied. The effect of the phase {(pA^ + (pp.) in 
the squark contributions is negligible since the sbottom mixing is small and moreover |Xfe| ^ 
\n\ tan/5, making its phase dependence insignificant. The slepton contributions (entering via 
the diagrams in Fig. 2) yield a shift in Mw of up to 10 MeV for small M?, i.e. about a 
third of the squark contributions. Even for the slepton contributions the dominant effect 
(about 60% of the total shift in Mw) can be associated with Ap, as a consequence of the 
D-term splitting of the sleptons. For large M? the slepton contributions show the expected 
decoupling behaviour. 
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Figure 10: Squark contributions to 6M]y as function of the phase (0^ + 



^M/' 



where (i>A = 



for different values of the common sfermion mass M 



other relevant SUSY parameters are set to tan/5 = 5, \Ath\ 



^ _ 500, 600, 1000 GeV. The 
2Mj, |/i| = 900 GeV. 



4.1.2 Chargino and neutralino sector dependence 

In this subsection we analyse the one-loop contributions from neutralinos and charginos to 
6Mw, entering via the self-energy, vertex and box diagrams shown in Figs. 6 and 7. In this 
sector the parameters Mi, M2 and /i can be complex. However, there are only two physical 
complex phases since one of the two phases of Mi and M2 can be rotated away. As commonly 
done we choose to rotate away the phase of M2. Generally, the phase dependence in this 
sector can be expected to be smaller than in the sfermion sector since the chargino/neutralino 
mass matrices are dominated by their diagonal elements, and the mass eigenvalues are mainly 
determined by |/x|, |Mi_2|, so that their phase dependence is small. 

The phase dependence is illustrated by Fig. 14, where the chargino/neutralino contribu- 
tions to Mw are shown as function of |yu| (left panel) and as contour lines in the 0^-|/i| plane 
(right panel) for 0Mi = 0. The other parameters are M? = 500 GeV, Mi = M2 = 200 GeV 
and tan/3 = 10. The effect of varying (p^ is much smaller than the overall contribution of the 
chargino/neutralino sector. In the scenario of Fig. 14 the chargino/neutralino contributions 
lead to a shift in the prediction for M^ of up to 11 MeV for small \fi\, while the effect of 
varying 0^ does not exceed 1 MeV. We have checked that also the dependence on (pMi is 
insignificant. 

In Fig. 15 we investigate the impact of varying \Mi\ and IM2I for zero complex phases. 
The other relevant parameters are M? = 250 GeV, |yu| = 300 GeV, tan/3 = 10. The shift in 
M-\v induced by varying IM2I can reach up to 15 MeV (i.e. the anticipated LHC precision). 
This is larger than the maximum shift in Fig. 14 because of the smaller sfermion masses. On 
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Figure 11: Contour lines of the squark contributions to 6Mw in the plane of 
|/i|, where (pA = 4>At = 4'Af The left plot shows a scenario with tan/? = 5, 
|A,b| = 1000 GeV, while in the right plot tan/3 = 30, Mj = 600 GeV, |A,6 



A + 0^) and 
Mf = 500 GeV, 
= 1200 GeV. 



the other hand, the effect of varying |Mi| stays below ~ 2 MeV. 

4.1.3 Gauge boson and Higgs sector dependence 

Finally we discuss the one-loop effects on Mw from the gauge and Higgs sectors. In order 
to identify the genuine SUSY effects, we compare the contribution in the MSSM with the 
one in the SM (where the Higgs-boson mass in the SM is set equal to the mass of the light 
CV-even Higgs boson of the MSSM, Mjjsm = Mjf^^^). The corresponding diagrams are 
shown in Figs. 3, 4 and 5, where only the diagrams shown in Fig. 4 differ in the MSSM and 
the SM. The parameters governing the Higgs sector at the tree level are Ma and tan/5 in 
the MSSM and M^sm in the SM. Since we incorporate higher-order contributions into the 
predictions for the MSSM Higgs masses and mixing angles, which are evaluated using the 
program FeynHiggs, further SUSY parameters enter the prediction for the gauge-boson and 
Higgs sector contribution. The effect of complex phases entering via the MSSM Higgs sector 
is formally of two-loop order. We therefore restrict to real parameters in this subsection. 

In Fig. 16 the shift 6Mw is given as a function of Ma in the MSSM and in the SM (with 
Mhsm = M^ss^'). The parameter tan/5 is fixed to tan/5 = 5,25, which affects M^^^^ 
and accordingly also Mjjsm. The other SUSY parameters are chosen as M? = 600 GeV, 
A,b = 1200 GeV, /i = 500 GeV, m~g = 500 GeV, Ms = 500 GeV. Fig. 16 shows that the 
overall effect of the gauge and Higgs boson sector is rather large, up to —60 MeV, which 
corresponds to about twice the current experimental error on Mw- The shift in the W boson 
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Figure 12: Squark contributions to 5 My/ as function of a common sfermion mass M? for 
different values of (0^* + 4'^l)■ The other relevant SUSY parameters are: \At^b\ = 350 GeV, 
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Figure 14: Neutralino and chargino contributions to 6Mw. The left plot shows 6M\y as 
function of |/i| for different values of (p^, while the right plot displays the contour lines of 
SM]v in the 4>^-\fi\ plane. The other relevant SUSY parameters are: (pMi = 4>M2 = 0, 
|Mi| = iMal = 200 GeV, tan (3 = 10 and M? = 500 GeV. 



mass obtained in the MSSM is slightly larger than in the SM. However, this genuine SUSY 
effect, i.e. the difference between the MSSM and the SM value of 6Mw, is always below 
~ 2 MeV for tan/9 = 5 and smaller than ~ 5 MeV for tan/3 = 25. The genuine SUSY effect 
is therefore below the anticipated GigaZ accuracy. 



4.2 Prediction for Mw 

In this section we combine all known contributions as described in Sect. 3.2 and analyse the 
parameter dependence of this currently most accurate MSSM prediction for M^ in various 
scenarios. These predictions are compared with the current experimental value for M^ [2], 



M^P 



80.404 ± 0.030 GeV 



(42) 



4.2.1 Dependence on SUSY parameters 



We start by comparing our full MSSM prediction for Mw with the corresponding SM value 
(with MffSM = Mh) as a function of Ma in Fig. 17. Like in Fig. 16 the other parameters 
are M^ = 600 GeV, A,6 = 1200 GeV, /i = 500 GeV, m~g = 500 GeV, Ms = 500 GeV. tan/3 
is set to tan/3 = 5,25. It can be seen in Fig. 17 that for this set of parameters the MSSM 
prediction is about 20 MeV higher than the SM prediction. While the MSSM prediction is 
within lo" of the experimental value of Mw, the SM prediction lies in the l-2o" interval. 
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Figure 15: Neutralino and chargino contributions to 5Mw as function of |Mi| and IM2I. The 
other relevant SUSY parameters are: 0a/^ = 0Af2 = 0m ~ 0' l/^l ~ ^^O GeV, tan/3 = 10 and 
Mj = 250 GeV. 



In Fig. 18 we show the prediction for M^/ as a function of M? and indicate how this 
prediction changes if the top-quark mass is varied within its experimental la interval, nit = 
(172.5 ± 2.3) GeV [58]. The other parameters are At,b = 2Mj, ^ = Ma = m-g = M2 = 
300 GeV, and tan/3 = 10. The result is compared with the current experimental value of 
M^. It can be seen in Fig. 18 that the observable M^ exhibits a slight preference for a 
relatively low SUSY scale, see also Refs. [5, 6] for a recent discussion of this issue. For the 
current experimental central value of the top-quark mass, Mw lies in the experimental la- 
interval only for a SUSY mass scale of M? ^ 800 GeV. Increasing rrit by one standard 
deviation allows M? up to at least 1300 GeV at the la level for this set of SUSY parameters. 

In Fig. 19 we show the dependence of M^ on /i and M2. The other parameters are set to 
Mj = 300 GeV, A,6 = 2 M^-, Ma = 1000 GeV, rUg = 600 GeV, tan/3 = 10. As can be seen in 
Fig. 19, varying /i between about 200 GeV and 1000 GeV results in a downward shift of more 
than 40 MeV in Mw- This strong dependence on the /i parameter is due to the neutralino 
and chargino as well as the squark contributions (one- and two-loop) to Ar. The sensitivity 
to the /i parameter from both MSSM particle sectors adds up and leads to the large shifts 
shown in Fig. 19. The neutralino and chargino contributions are responsible for roughly one 
third of the 40 MeV shift for small M2 and become negligible for large M2, with the remaining 
MSSM parameters specified as given above. The shift induced by varying M2 between about 
100 GeV and 1000 GeV (as explained above, M2 and Mi are varied simultaneously according 
to eq. (37)) amounts up to about 15 MeV in M^y. For the relatively small value chosen for 
the common sfermion mass, M? = 300 GeV, all combinations of fi and M2 yield a result 
within la of the experimental result of M-^. Using instead a larger SUSY mass scale of, for 
instance, M? = 600 GeV would result in Mw values within the experimental la interval only 
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Figure 16: The contribution 6Mw to the prediction for the W boson mass from the gauge- 
boson and Higgs sector contributions in the MSSM and the SM, where the SM contributions 
are evaluated for Mjjsm = Mjf^^^. The shift 6M]y is given as a function of the CP-odd 
Higgs-boson mass M^ and tan/? = 5,25. The other SUSY parameters were chosen to be: 

1200 GeV, /i = 500 GeV, m. = 500 GeV, M2 = 500 GeV. 



Mj = 600 GeV, At,b 



for /i < 500 GeV. 

Finally we discuss the effect of varying the complex phase (pAt on the prediction for Mw- 
As explained above, we use our complete one-loop result for the phase dependence and employ 
eq. (32) to approximate the effect of the complex phases at the two-loop level. In Fig. 20 the 
prediction for M^ is shown as a function of (pAt for |A,b| = 1000 GeV, (pA^ = 0, tan/? = 10, 
m-g = 500 GeV, Ma = 500 GeV, M2 = 250 GeV and /i = 500 GeV. The results are plotted 



The dependence on (pAt is at most of the order 
For heavier sfermions with Mj = 1000 GeV only a 



for Mf = 500 GeV, 600 GeV and 1000 GeV. 
2 MeV for Mj = 500 GeV, 600 GeV. 
1 MeV shift in Mw can be observed. 

As explained above, our result for Mw goes beyond the results previously known in the 
literature [3] because of the inclusion of complex phases and an improved treatment of higher- 
order SM contributions. We have checked that for real parameters our new result agrees with 
the previously most advanced implementation [3] typically within about 5 MeV. 



4.2.2 The SPS inspired benchmark scenarios 

In this subsection we show Mw in the SPS la, SPS lb and SPS 5 benchmark scenarios [34]. 
This should give an indication of the Mw prediction within "typical" constrained MSSM 
(CMSSM) scenarios. In the original definition the SPS parameters are DR parameters. Here 
we treat them as on-shell input parameters for simplicity, since the effects of the DR to on- 
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Figure 17: Prediction for Myy in the MSSM and the SM, where the SM contributions are 
evaluated for M^sm = Mjf^^^. The prediction for Mw is shown as function of the CP-odd 
Higgs-boson mass M^ for tan/5 = 5,25. The other SUSY parameters are: M? = 600 GeV, 
At,b = 2 M^-, /i = 500 GeV, m-g = 500 GeV, Ma = 500 GeV. 
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Figure 18: Prediction for M^ as function of a common sfermion mass for nit = (172.5 ± 
2.3) GeV. The other SUSY parameters are set to be: A^ ^ = 2 M?, tan/3 = 10, /x = 300 GeV, 
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Figure 19: Prediction for Mw as function of M2 for /i = 250 GeV, 500 GeV, 1000 GeV (left 
plot) and as function of /i for M2 = 150 GeV, 250 GeV, 500 GeV, 1000 GeV (right plot). 



The other SUSY parameters are: M? = 300 GeV, At^b 
Ma = 1000 GeV. 
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Figure 20: Prediction for Mw as a function of the phase of At using relation (32). The SUSY 
parameters are \At,b\ = 1000 GeV, ^a^ = 0, tan^ = 10, m-g = 500 GeV, Ma = 500 GeV, 
M2 = 250 GeV and /i = 500 GeV. 



shell transition are expected to be small and therefore irrelevant for the qualitative features 
that we discuss. In order to analyse the dependence of M\y on the scale of supersymmetry 
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we introduce a scale factor. Every SUSY parameter of mass dimension of the considered SPS 
point is multiplied by this parameter, i.e. Ma = (scalefactor) x M^^, Mp p, = (scalefactor) x 
M?^?,, AtM = (scalefactor) x A^J^, fi = (scalefactor) x fi^^^, Mi^2,3 = (scalefactor) x Mf g,!- 
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Figure 21: Prediction for Mw within the SPS la, SPS lb and SPS 5 inspired scenarios. 
Mw is shown as a function of mi^, the lighter of the two stop squarks (left plot), and as a 
function of Ma, the mass of the CP-odd Higgs boson (right plot). The SPS parameters of 
mass dimension are varied with the scale of supersymmetry as described in the text. 



In Fig. 21 we show the result for the three SPS scenarios as a function of the lighter t 
mass, rrii^^ (left plot), and as a function of Ma (right plot). The prediction for Mw is similar 
in all three scenarios. The variation of the dimensionful SUSY parameters shifts the MSSM 
prediction for Mw by up to 35 MeV. As can be seen in the left plot of Fig. 21, agreement 
at the la level with the experimental result is obtained for m^^ < 600 GeV. Since we scale 
all dimensionful parameters simultaneously, the variation from small to large m^^ (left plot) 
is the same as the one from small to large Ma (right plot). However, for the same Ma value 
the three scenarios can differ also by up to ~ 30 MeV. 

For large values of the SUSY mass scale one expects a decoupling behaviour of the SUSY 
contributions, i.e. one expects that the prediction for Mw in the MSSM coincides with the SM 
prediction (for Mji^sm = M^^^^) in the limit of large SUSY masses. We analyse the decou- 
pling behaviour in Fig. 22 for the SPS la scenario. We compare the MSSM prediction with 
the corresponding SM prediction of Mw with Mjjsm = Mjf-^^^. For mj:_^ ^ 500 GeV large 
deviations between the MSSM and the SM prediction can be observed. For m^^ ^ 2500 GeV 
the difference drops below the level of 1 MeV, and for even larger m^^ values the MSSM 
result converges to the SM result. It should be noted in this context that the prescription 
described in Sect. 3.2.1 (see eq. (29)) has been crucial in order to recover the most up-to-date 
SM prediction for Mw in the decoupling limit. 
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Figure 22: Prediction for M^ in the decoupling limit compared to the SM prediction with 
MfjSM = Mjf^^^. My/ is shown as function of the light stop mass mi^ in the SPS la 
benchmark scenario. All dimensionful SUSY parameters are scaled with a common factor, 
as described in the text. 



4.2.3 The focus point scenario 

As we have seen in the previous section, a relatively light SUSY scale leads to a prediction 
for Mw within the CMSSM (and of course also the unconstrained MSSM) that is in slightly 
better agreement with the experimental value of M\y than the SM prediction. A region of the 
CMSSM that has found a considerable interest in the last years is the so-called focus point 
region [67]. This region is characterized by a relatively small fermionic mass parameter, r7ii/2, 
while the common scalar mass parameter itiq is very large, and also tan/3 is relatively large. ^ 
We now investigate whether it is also possible to obtain a prediction for Mw that is in better 
agreement with the experimental value than the SM prediction if the MSSM parameters are 
restricted to the focus point region. 

We have evaluated three representative scenarios, using rrit = 172.5 GeV, tan/3 = 50 and 
/x > 0. We have chosen a point with the currently lowest value of mi/2 in the focus point 
region for which the dark matter density is allowed by WMAP and other cosmological data 
(see e.g. Ref. [6] for a more detailed discussion) and two further points with higher mi/2 
(and higher mo) along the strip in the 'mi/2-mQ plane that is allowed by the dark matter 
constraints. These points yield the following results for M^ in the MSSM 



"^1/2 = 

^ Mw 



250 GeV, mo = 1650 GeV, Aq = -250 GeV 



80.380 GeV, M^{M, 



//SM 



M, 



MSSM\ 
h ) 



80.361 GeV, 



(43) 



^The CMSSM is characterized in terms of three GUT-scale parameters, the common fermionic mass 
parameter mi/2, the common scalar mass parameter mo, and the common trihnear couphng Aq. These 
high-scale parameters are supplemented by the low-scale parameter tan/3 and the sign of /i. 
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mi/2 = 330 GeV, rriQ = 


2030 GeV, A^ = 


-330 GeV 






^Mw = 80.372 GeV, 


MI^{Mhsm 


= M.^SSM) ^ 


: 80.360 GeV, 


(44) 


mi/2 = 800 GeV, rriQ = 


3685 GeV, Aq = 


-800 GeV 






^ Mw = 80.361 GeV, 


MI^{Mhsm 


= Mf ss^i) = 


: 80.359 GeV, 


(45) 



(2) 
(3) 



where the low-scale parameters of the MSSM have been obtained from the high-scale pa- 
rameters mi/2, mo, Ao with the help of the program ISAJET 7.71 [68]. For comparison, the 
corresponding prediction in the SM with Mjjsm = M^^^^ is also given. 

One can see from eqs. (43)-(45) that only for the point with the lowest 1711/2 value a 
large difference of up to ~ 20 MeV occurs in comparison to the SM result. Such low mi/2 
values correspond to very low masses of neutralinos and charginos. As a consequence, the 
main contribution to the shift in Mw arises from the chargino and neutralino sector (see 
Sect. 4.1.2). Since for the ttiq value of eq. (43) the squarks are not completely decoupled, 
there is also a small contribution from the squark sector. Already with slightly higher mi/2 
values, eq. (44), the contribution to Mw becomes much smaller. For even higher values, 
eq. (45), the resulting prediction for Mw is very close to the corresponding SM (decoupling) 
limit. These predictions deviate by about 1.5o" from the current experimental value of Mw, 
eq. (42). The focus point region therefore improves the prediction for Mw only for very low 
mi/2- For most of the allowed parameter space, however, the improvement is small. The 
deviation in Mw contributes to the relatively bad fit quality of the focus point region in a fit 
to electroweak precision observables as observed in Ref. [6]. 

4.2.4 Split SUSY 

Another scenario that has recently found attention is the so-called "split-SUSY" scenario [69] . 
Here scalar mass parameters are made very heavy and only the fermionic masses (i.e. the 
chargino, neutralino and gluino masses) are relatively small. According to the analysis in 
the previous sections only a small deviation in the Mw prediction from the SM limit is to be 
expected. 

We have evaluated the prediction for Mw in the split-SUSY scenario. In Fig. 23 we show 
the SUSY contribution to Mw, i-e. the deviation between the MSSM result in the split-SUSY 
scenario and the corresponding SM result. This deviation is obtained by choosing a large 
value for M? and subtracting the SM result with Mj^sm = M^^^^. For definiteness we 
have chosen M? = 3 TeV and Ma = 2 TeV. Choosing a higher scale for the scalar mass 
parameters would lead to even slightly smaller deviations from the SM prediction than the 
ones shown in Fig. 23. The resulting shifts in Mw are displayed in Fig. 23 in the /X-M2 plane 
for tan/3 = 10. The gluino mass has been set to rrig = 300 GeV, however no visible change in 
Fig. 23 occurs even for rrig = 3000 GeV. As expected, only for rather small values of /i and 
M2, M2 < 400 GeV and |/i| < 500 GeV a deviation from the SM limit larger than 5 MeV 
is found (a similar result has been obtained in Ref. [70], see also Ref. [71]). Even with the 
GigaZ precision for Mw only a very light chargino/neutralino spectrum would result in a lo" 
deviation in Mw compared to the SM prediction. 
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Figure 23: Difference between the result for M^J\/ in tlie MSSM and the SM for large scalar 
fermion masses. The SUSY parameters are M? = 3000 GeV, A^ ^ = 2 M ?, rrig = 300 GeV, 
M^ = 2000 GeV and tan/3 = 10. 



4.2.5 MSSM parameter scans 

Finally, we analyse the overall behaviour of Mw in the MSSM by scanning over a broad range 
of the SUSY parameter space. The following SUSY parameters are varied independently of 
each other, within the given range, in a random parameter scan: 



sleptons 


Mpp, = 100 . . . 2000 GeV 


light squarks 


Mpp, = 100 ... 2000 GeV 

' up/down 


t/b doublet 


Mpp, = 100 ... 2000 GeV 

' up/down 




At^h = -2000 . . . 2000 GeV 


gauginos 


Mi,2 = 100 . . . 2000 GeV 



Higgs 



mg = 195... 1500 GeV 
/i = -2000 . . . 2000 GeV 
Ma = 90 - 1000 GeV 
tan/5 = 1.1... 60 



(46) 



We have taken into account the constraints on the MSSM parameter space from the LEP 
Higgs searches [72,73] and the lower bounds on the SUSY particle masses from Ref. [35]. 
Apart from these constraints no other restrictions on the MSSM parameter space were made. 
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In Fig. 24 we show the result for M^ as a function of the hghtest t mass, mf:_^. The 
top-quark mass has been fixed to its current experimental central value, rrit = 172.5 GeV. 
The results are divided into a dark (green) shaded and a light (green) shaded area. In the 
latter at least one of the ratios rn^Jm^^ or rni^/mi^ exceeds 2.5,"^ i.e. the darker shaded region 
corresponds to a moderate splitting among the t or 5 doublets. In this region the MSSM 
prediction for Mw does not exceed values of about 80.550 GeV. In the case of very large 
splitting in the t and h doublets, much larger My/ values up to 81.150 GeV would be possible 
(which are of course ruled out by the experimental measurement of M\y). 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



81.15 




M „ "^ = 80.404 + 0.030 GeV 



250 500 750 1000 1250 1500 1750 

m~ [GeV] 

Figure 24: Prediction for My/ as function of m^^, the mass of the lighter stop squark. The 
SUSY parameters are varied independently of each other in a random parameter scan as 
described in the text. The top-quark mass is fixed at its current experimental central value, 
mt = 172.5 GeV. 



In Fig. 25 we compare the SM and the MSSM predictions for Myy as a function of mt as 
obtained from the scatter data. The predictions within the two models give rise to two bands 
in the rrit-Myr plane with only a relatively small overlap region (indicated by a dark-shaded 
(blue) area in Fig. 25). The allowed parameter region in the SM (the medium-shaded (red) 
and dark-shaded (blue) bands) arises from varying the only free parameter of the model, 
the mass of the SM Higgs boson, from MfjsM = 114 GeV, the LEP exclusion bound [73] 
(upper edge of the dark-shaded (blue) area), to 400 GeV (lower edge of the medium-shaded 
(red) area). The very light-shaded (green), the light shaded (green) and the dark-shaded 

^ We work in the convention that rrif < mi . 

Jl — 12 
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(blue) areas indicate allowed regions for the unconstrained MSSM. In the very light-shaded 
region (see Fig. 24) at least one of the ratios rn^^/m^^ or Tn^^/mi_^ exceeds 2.5, while the 
decoupling limit with SUSY masses of 0{2 TeV) yields the lower edge of the dark-shaded 
(blue) area. Thus, the overlap region between the predictions of the two models corresponds 
in the SM to the region where the Higgs boson is light, i.e. in the MSSM allowed region 
(M/i < 135 GeV [38,39]). In the MSSM it corresponds to the case where all superpartners are 
heavy, i.e. the decoupling region of the MSSM. The current 68% CL. experimental results^ for 
rrit and Mw are indicated in the plot. As can be seen from Fig. 25, the current experimental 
68% CL. region for rrit and My/ exhibit a slight preference of the MSSM over the SM. The 

,• • f iU TH , /T TTn /r Tcvatron/LHC i ri t r r 71 j-Tcvatron/LHC 

prospective accuracies tor the ievatron/LML [pm^ = 1 GeV, oM^r = 

15 MeV) and the ILC with GigaZ option ((5mJ^^/^^^^^ = 0.1 GeV, 6M^^^^^'^^^ = 7 MeV) 
are also shown in the plot (using the current central values), indicating the potential for a 
significant improvement of the sensitivity of the electroweak precision tests [74] . 

4.3 Remaining higher-order uncertainties 

As explained above, we have incorporated all known SM corrections into the prediction for 
Mw in the MSSM, see eq. (29). This implies that the theoretical uncertainties from unknown 
higher-order corrections reduce to those in the SM in the decoupling limit. In the SM, based 
on all higher-order contributions that are currently known, the remaining uncertainty in Myy 
has been estimated to be [16] 

SM^ = 4 MeV . (47) 

Below the decoupling limit an additional theoretical uncertainty arises from higher-order cor- 
rections involving supersymmetric particles in the loops. This uncertainty has been estimated 
in Ref. [30] for the MSSM with real parameters depending on the overall sfermion mass scale 

SMw = 8.5 MeV for Mj < 500 GeV, 

SMw = 2.7 MeV for Mj = 500 GeV, (48) 



SMw = 2.4 MeV for M,- = 1000 GeV. 



The full theoretical uncertainty from unknown higher-order corrections in the MSSM with 
real parameters can be obtained by adding in quadrature the SM uncertainties from eq. (47) 
and the SUSY uncertainties from eq. (48). This yields 6Mw = (4.7 — 9.4) MeV depending 
on the SUSY mass scale [30]. 

Allowing SUSY parameters to be complex adds an additional theoretical uncertainty from 
unknown higher-order corrections to the Mw prediction. While at the one-loop level the full 
complex phase dependence is included in our evaluation, it is only approximately taken into 
account at the two-loop level as an interpolation between the known results for the phases 
and n, see eq. (32). We concentrate here on the complex phase in the scalar top sector, (pAt 
(we keep 0^ fixed) , since our analysis above has revealed that the impact of the other phases 
is very small already at the one-loop level. 

We estimate the uncertainty from unknown higher-order corrections associated with the 
phase dependence as follows. The full result for Mw{4>), i-e. for a given phase (where the 
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The plot shown here is an update of Refs. [3,24]. 
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Figure 25: Prediction for Miy in the MSSM and the SM as a function of rrit in comparison 
with the present experimental results for M^ and rrit and the prospective accuracies (using 
the current central values) at the Tevatron / LHC and at the ILC The allowed region in 
the MSSM, corresponding to the light-shaded (green) and dark-shaded (blue) bands, results 
from varying the SUSY parameters independently of each other in a random parameter scan. 
The allowed region in the SM, corresponding to the medium-shaded (red) and dark-shaded 
(blue) bands, results from varying the mass of the SM Higgs boson from Mh = 114 GeV to 
Mh = 400 GeV. Values in the very light shaded region can only be obtained in the MSSM 
if at least one of the ratios mi^/mi^ or mi^/mi^ exceeds 2.5. 



two-loop corrections are taken into account as described in Sect. 3.2.3) lies by construction 
in the interval 

[Mr(0) + (M^^(0) - M^L(O)), Mr(vr) + {M'^{<P) - M^L(7r))] . (49) 

The minimum difference of M^^'((/)) to the boundary of this interval, 

Min{|Mr(0)-(Mr(O) + (M^L(0)-M^L(O)))|,|Mr(0)-(AC(7r) + (M^L(0)-M^L(vr)))|}, 

(50) 
can be taken as estimate for the theoretical uncertainty (this automatically ensures that no 
additional uncertainties arise for cj) = 0, tt). 

As representative SUSY scenarios we have chosen SPS la, SPS lb, and SPS 5, each for 
Mj = 1000 GeV, 500 GeV, and for Mj < 500 GeV ^ (as above we have varied Mj, At,b and 

^The lowest values considered for Mf are roughly 300, 300, 400 GeV for SPS la, SPS lb, SPS 5, 
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5Mw -- 


= 3.2 MeV for Mj 


< 500 GeV, 


SMw = 


= 2.0 MeV for Mj 


= 500 GeV, 


SMw -- 


= 0.7 MeV for Mj 


= 1000 GeV. 



IJ, using a common scale factor). A non-zero complex phase (pAt has been introduced as an 
additional parameter. In order to arrive at a conservative estimate for the intrinsic error, we 
use the value obtained for (p^t = ±7r/2. In this way we obtain 



(51) 



The full theoretical uncertainty from unknown higher-order corrections in the MSSM with 
complex parameters can now be obtained by adding in quadrature the SM uncertainties from 
eq. (47), the theory uncertainties from eq. (48) and the additional SUSY uncertainties from 
eq. (51). This yields 6Mw = (4.7 — 9.9) MeV depending on the SUSY mass scale. 

The other source of theoretical uncertainties besides the one from unknown higher-order 
corrections is the parametric uncertainty induced by the experimental errors of the input 
parameters. The current experimental error of the top-quark mass [58] induces the following 
parametric uncertainty in M^ 

(^mfP = 2.3 GeV => SM^'^^""' = 14 MeV , (52) 

while the uncertainty in Aa^^^i [75] results in 

5(Aa[2j) = 3.6 X 10-^ ^ 5M^""''^°''^^<^ = 6.5 MeV . (53) 

The uncertainty in irit will decrease during the next years as a consequence of a further 
improvement of the accuracies at the Tevatron and the LHC. Ultimately it will be reduced by 
more than an order of magnitude at the ILC [76]. For AaJ^j^^ one can hope for an improvement 
down to 5 X 10^^ [77], reducing the parametric uncertainty to the 1 MeV level (for a discussion 
of the parametric uncertainties induced by the other SM input parameters see e.g. Ref. [3]). 
The effect of 5{/S.a\^l^ on Fig. 25 is small. In order to reduce the theoretical uncertainties from 
unknown higher-order corrections to the 1 MeV level, further results on SM-type corrections 
beyond two-loop order and higher-order corrections involving supersymmetric particles will 
be necessary. 

5 Conclusions 

We have presented the currently most accurate evaluation of the W boson mass in the 
MSSM. The calculation includes the complete one-loop result, taking into account for the 
first time the full complex phase dependence, and all known higher-order corrections in the 
MSSM. Since the evaluation of higher-order contributions in the SM is more advanced than 
in the MSSM, we have incorporated all available SM corrections which go beyond the results 
obtained so far in the MSSM. Our prediction for My/ in the MSSM therefore reproduces the 
currently most up-to-date SM prediction for My/ in the decoupling limit where the masses of 



respectively. For lower values the parameter points are excluded by Higgs mass constraints. The light stop 
mass for the SPS 5 point lies considerably below 400 GeV. 
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all supersymmetric particles are large. A public computer code based on our result for M]y 
is in preparation. 

We have analysed in detail the impact of the various sectors of the MSSM on the prediction 
for M\Y, focussing in particular on the dependence on the complex phases entering at the one- 
loop level. The most pronounced phase dependence occurs in the stop sector, where the effect 
of varying the complex phase that enters the off-diagonal element in the stop mass matrix can 
amount to a shift of more than 20 MeV in Mw- It should be noted, however, that the complex 
phases in the squark sector at the one-loop level enter only via modifications of the squark 
masses and mixing angles. As a consequence, a precision measurement of the CP-conserving 
observable M^y alone will not be sufficient to reveal the presence of CP-violating complex 
phases. The phase dependence of M^ will however be very valuable for constraining the 
SUSY parameter space in global fits where all accessible experimental information is taken 
into account. 

We have illustrated the sensitivity of the precision observable Mw to indirect effects of 
new physics by comparing our MSSM prediction with the SM case. Confronting the MSSM 
prediction in different SUSY scenarios and the SM result as a function of the Higgs-boson 
mass with the current experimental values of Mw and rrit, we find a slight preference for 
non-zero SUSY contributions. As representative SUSY scenarios we have studied various SPS 
benchmark points, where we have varied all mass parameters using a common scalefactor. 
The MSSM prediction lies within the la region of the experimental Mw value if the SUSY 
mass scale is relatively light, i.e. lower than about 600 GeV. 

The prospective improvement in the experimental accuracy of Mw at the next generation 
of colliders will further enhance the sensitivity to loop contributions of new physics. We 
have found that even for a SUSY mass scale of several hundred GeV the loop contribution 
of supersymmetric particles to Mw is still about 10 MeV, which may be probed in the high- 
precision environment of the ILC In SUSY scenarios with even higher mass scales, however, 
it is unlikely that SUSY loop contributions to Mw can be resolved with the currently foreseen 
future experimental accuracies. We have studied in this context the focus point and the split- 
SUSY scenarios, which have recently received significant attention in the literature. In the 
focus point scenario only at the lower edge of the allowed parameter region in the mi/2-mo 
plane a sizeable contribution to Mw can be achieved. For higher values, as well as for the 
split-SUSY scenario, we find that the SUSY loop effects are in general very small, so that it 
is not possible to bring the prediction for Mw significantly closer to the current experimental 
central value as compared to the SM case. 

Finally we have analysed the theoretical uncertainties in the Mw prediction that arise 
from the incomplete inclusion of the complex phases at the two-loop level. We estimate that 
this uncertainty can amount up to roughly 3 MeV in the prediction for Mw, depending on the 
SUSY mass scale. Combined with the estimate of the possible effects of other unknown higher- 
order corrections we find that the total uncertainty from unknown higher-order contributions 
is currently about 10 MeV for small SUSY mass scales. 
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